The effect of material parameter values on 
the relation between the energy gap width 
and the scatterer symmetry in 
two-dimensional phononic crystals 

Izabela Sliwaf, Maciej Krawczyk| 

f Nonlinear Optics Division, Adam Mickiewicz University, Umultowska 85, 
61-614 Poznan, Poland 

\ Surface Physics Division, Adam Mickiewicz University, Umultowska 85, 61-614 
Poznan, Poland 

E-mail: krawczyk9amu.edu.pl 

Abstract. The phononic band structures of two-dimensional solid phononic 
crystals with different lattice and scatterer symmetry are studied numerically, 
with three types of lattice (square, triangular and rectangular) and four different 
scatterer shapes (circle, hexagon, square and rectangle) considered. XY and 
Z vibration modes are investigated separately. Two types of phononic crystal 
are considered: one composed of high-density rods embedded in a low-density 
matrix, the other of low-density rods in a high-density matrix. In the former 
case, lattice type and polarization being fixed, the broadest gaps are obtained 
when the symmetry of the rods corresponds to that of the lattice (the shape of 
a rod is identical with that of the first Brillouin zone); the largest gap width 
values are observed in triangular lattice-based crystals (compared to those based 
on square and rectangular lattices), the shape of the corresponding first Brillouin 
zone being closest to a circle. These rules do not apply to structures in which the 
density of the rod material is lower than that of the matrix. In this case, when 
the symmetry of the rods corresponds to that of the lattice, gaps either fail to 
appear at all, or are much narrower than in other configurations. The effect of 
other material parameter values (such as the longitudinal and transversal velocity 
values) on the considered relation proves much lesser. 
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1. Introduction 

Elastic wave propagation in periodic structures referred to as phononic crystals (PC) 
has been the subject of intensive research over the last few years P El EI El El Cj • 
Composed of elastic centers (e. g. cylinders, in the two-dimensional case) embedded 
in a matrix of a different elastic material and disposed in nodes of a crystal lattice, 
phononic crystals show in their band structure energy gaps forbidden to elastic 
wave propagation. Due to this property, these structures have extensive practical 
applications, especially as so-called acoustic screens, suppressing waves from a given 
frequency range JUj. For this reason, it is important to know how to design structures 
with the phononic gap as wide as possible. 

Beside the physical properties of the PC crystal component materials, the 
phononic gap width is found to depend strongly on the lattice symmetry, as well 
as on the scatterer shape ^] El ■ This dependence is a common feature of phononic 
and photonic crystals p^l- In both crystals, the gap width depends also on the rod 
orientation with respect to the lattice, and thus on the relative orientation of the rods 

as weii Ha mj. 

The effect of scatterer shape and symmetry on the phononic band gap in two- 
dimensional solid phononic crystals has recently been studied by W.Kuang et al. 
[TTj : the investigated structures consisted of steel rods embedded in epoxy, and 
thus represented the case of high-density rods in a low-density matrix only. Three 
types of lattice (triangular, hexagonal and square) and four different scatterer shapes 
(hexagon, circle, square and triangle) were considered. For a given lattice symmetry, 
the absolute phononic band gap width was found to be the largest when the shape 
and orientation of the scatterers corresponded to those of the lattice. Besides, for 
a given scatterer shape, the gap was found to be the broadest when the lattice 
had the largest coordination number, as the crystal symmetries in this case were 
not reduced by the scatterers. Also, the band gap was shown to be controllable by 
adjusting the scatterer orientation and size: the normalized gap width was found to 
increase as a result of scatterer rotation if the latter caused the maximum allowed 
filling fraction to increase; otherwise the normalized gap width would be reduced. 
The aim of this study is to show the effect of rod and matrix material parameters 
(the mass density, in particular) on the relation between the energy gap width and 
the symmetry of the scatterers and the crystal lattice in two-dimensional phononic 
crystals. This objective will be realized by investigating the dependence of the gap 
width on the rod symmetry, with the rods embedded in different sorts of lattice; the 
two cases: high-density rods in low-density matrix, and low-density rods in high- 
density matrix, will be considered separately. Four composite types will be examined 
first: steel/epoxy, C/epoxy, Pb/epoxy and duralumin/epoxy, the epoxy component 
being always the matrix. Inverse structures, i.e. epoxy/steel, epoxy/C, epoxy/Pb 
and epoxy/duralumin, representing the case of low-density rods in a high-density 
matrix, will be investigated in the next step. Phononic band gaps will be studied in 
two-dimensional solid phononic crystals differing in lattice symmetry and rod shape. 
Specifically, three types of lattice: square (sq), triangular (tr) and rectangular (rc), and 
four different scatterer shapes: circle (C), hexagon (HX), square (SQ) and rectangle 
(RC), will be considered. The rod orientation with respect to the lattice is fixed, the 
symmetry axes of the rods corresponding to those of the crystal lattice with the same 
symmetry, the configuration allowing to generate the broadest energy gaps in phononic 
crystals composed of two solids |17| . In our investigation of elastic wave propagation 
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in these structures, XY modes (polarized in the wave propagation plane) and Z modes 
(polarized in the direction perpendicular to the wave propagation) will be examined 
separately. The effect of rod and matrix material parameters, especially the density 
ratio between the rod material and the matrix material, on the relation between the 
energy gap width and the lattice and rod symmetry proves to be substantial for both 
Z and XY modes. 

The paper is divided into four sections. In Section 2, we describe the plane wave 
(PW) method used in the band structure calculation; the numerical results, as well as 
a table of material parameters, are presented in Section 3; the results are discussed in 
Section 4. 



2. Method of calculation 



The phononic band structure is calculated through the elastic medium wave equation, 
solved with the use of the plane wave method ^1 E|. The plane wave method 
consists in transforming a differential wave equation, whose coefficients are periodic 
functions of two-dimensional position vector r = (x\,X2), into the reciprocal space. 
For a nonhomogeneous elastic medium the wave equation reads |lf>| : 



p(r) 



d 2 Uj 
dt 2 



V- [p{r)c 2 {f)Vui} +V 



p{r)c t (f) 



du 

dxi 



(1) 



d 



where Ct and c; are the elastic wave transversal and longitudinal velocities, respectively; 
p denotes mass density, and u(f, t) stands for the time and space-dependent deflection 
vector, describing the medium molecules deflection from the balance point, and defined 
by its three components it, (i = 1, 2, 3) in the Cartesian coordinate system (OX1X2X3). 

The structure is assumed to be homogeneous along the X3 axis, which reduces the 
problem to two dimensions. If the wave propagation is confined to the phononic crystal 
periodicity plane (the vector — (xi,X2) plane), splits into two independent 
equations [T5] : 

• an equation describing XY modes, or oscillations polarized in plane (X1OX2): 



p{r\\) 



If 



p(i)c?(i) 



dxi 



+ i [(/>(i) c '(i) - 2 <°(i) c ?(i)) Vt • u T ] 



(2) 



where 



Ut = u\e\ + u 2 e 2 , 

9 d _ 

Vt = q — ei + - — e2 
axi 0x2 

for i = 1, 2. This equation can be re-written with Lame coefficients: 



where 



and 



d 2 Ui 



/' 
A 



d 

dxi 



du 1 
dxi 



_d_ 

dxi 



/ duj 
\dxi 



dxi 



(3) 



pel 
P(c 2 



2c?) 
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• an equation describing Z modes, or oscillations polarized along the rod axis, i. e. 
along vector el (el _L k): 

P(n)^ = V-[p{r^{ri)Vu z ]. (4) 

As the coefficients in equations of motion <(3J) and 10} are periodic functions, 
through Bloch's theorem, the deflection vector, u(f, t), can be expressed as |12j : 

tf(rjj , t) = e^i £ %(G)e^"1 , (5) 

G 

k being a two-dimensional Bloch vector, and G denoting a reciprocal lattice vector. 
The reciprocal lattice vectors, G, for the lattices considered in this paper (Figure ^) 
are defined as follows: 



for a square lattice 



2tt 

G = del + del = — (niel + r^ei) , (6) 



• for a triangular lattice 

G = del + del 

• for a rectangular lattice 



27T — 

mv3ei + (— m + 2n 2 ) el 



a 



V3 



(7) 



2?r 

G = del + del = \n\Lie\ + n 2 Lie" 2 ] , (8) 

ni and n 2 being integers, a being a lattice constant in the case of square and triangular 
lattices, and L\ and L 2 being the rectangular lattice constants along the x\ and x 2 
directions, respectively. 

Those of the material parameters which are periodic functions can be Fourier- 
expanded. For Z modes, described by Q), Fourier expansion can be applied to the 
inverse of mass density, p~ x , and to coefficient /i; in @, describing XY modes, the 
Fourier-expandable parameters are \x and A: 

p- 1 (i) = E'" 1 ( < VH (9) 

G 

Kr|) = ^MG>^, (10) 

G 

A(rl)=^A(G)e^l (11) 

G 

Coefficients p~ 1 (G), p(G) and A(G) in the above expansions are calculated from 
the inverse Fourier transformation. All these coefficients will have the same form; for 
example, the inverse of mass density will read: 



P~ l (G) = - J J p-^r^e-^drl (12) 

the integration covering the two-dimensional unit cell surface, S. When G = 0, (11211 
defines the mean value: 

p-\G^0)^p^^p- A 1 f + p B 1 {l~f) 1 (13) 
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where pa is the mass density value in the rod, ps is its matrix counterpart, and 
/ = S r /S denotes the filling fraction (S r being the rod cross-section area). When 
G^O, {njl becomes: 

p-\G * 0) = {p~ A l - pz 1 ) 1(G) = Ap-'liG), (14) 
where Ap _1 = p^ 1 — Pg 1 , and 1(G) is a structural factor determined by the rod shape: 



W-8 



-iG-r, 



dflt 



(15) 



The integration in (|15fl is performed over the rod cross section. For the four types 
of rods (Figure with circular, hexagonal, rectangular or square cross section, the 
integration results in the following formulae (Ic(G), Ihx(G), Irc(G) and Isq(G) 
denoting the respective structural factors): 

• for cylindrical rods: 

,Ji(Gr ) 



Ic(G) = 2f- 



Gr 



(16) 



where J\ is a first-order Bessel function, r is the cylinder cross-section radius, 
and filling fraction / is defined as follows: 

2 

/ = — < / < — for a square lattice, 
a 1 4 

/ = j=, < / < —= for a triangular lattice, 

a 2 V3 V3 



and 



/ = 

for hexagonal rods: 



L1L2 



for a rectangular lattice; 



(G) = 



2/ 

3d; 



3d; 
/ 

G 2 l 



• l/^i 7\ i 1 — cos(dO 

sm(GiZ) + 



, d ^ 0, G 2 = 0, 



l-cos(2di) 



, sin(2Gi0 , |Gx| = |G 2 |/V3, G 2 ^ 0, 



cos(G2;/v / 3-d;)-cos(2G 2 ;/v / 3) 
(G 2 +d \/3); 

+ |G 2 |/V3, |G 2 |^0, 



S (G 2 ;/v / 3+Gi;)-cos( 2 G 2 ;/v / 3) 

(g 2 -GiV3); 



(17) 



where / = &-\/3/2, & denoting the regular hexagon side length; the corresponding 
filling fraction is: 

/=^5 0</<^ for a square lattice, 

36 2 

/ = — — , < / < 1 for a triangular lattice, 



and 



35V3 
2£iL 2 



for a rectangular lattice; 
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• for rectangular rods: 

f Isiii^), Gi = 0, G 2 ^0, 
^c(G) = ^ |f sin(^), G 2 = 0, Gi ^ 0, 

I g^ 8 ^^ 1 ) 8 ^^), Gx^O, G 2 ^0, 
where Zi and Z 2 are the rectangle side lengths, and filling fraction / reads: 

for a square lattice, 



(18) 



f= hh 
1 a 2 



and 



/ = 



/ 



2hh_ 
a 2 V3 

hh 
L1L2 



for a triangular lattice, 



for a rectangular lattice; 



• the structural factor for square rods, Isq(G), is deduced from lfT%|l through 
putting h = h = I' 



Isq(G)=I rc (G) with h = h = l. 



(19) 



The maximum filling fraction values for rectangular lattice and rods depend on the 
rectangle side length ratio, l x /l y and the lattice constant ratio, L x /L y . 

The substitution of 10 -EJ, CBt and (fT4*|) into J3J) leads to an infinite system 
of algebraic equations for XY mode eigenvalues ui(k): 



u 2 u l r -. 

k+G 



E G - {Elg<> P^(G - G") [KG" ~ G')(k + G")i(k + G")i 
+/x(G" - G')(k + G')i(k + G")i 



u Ug> 



(20) 



- £ G ~„ [p- 1 (G - G>')m(G>' - G') Ej (k + G')j (k + G'%] ul +& } , 



where i,j, I — 1,2. Analogically, the following Z-mode eigenproblem is deduced from 
Q: 



w 2 u E (G')= P^iG-G")^" -G')(k + G')- (k + G")u % {G'). 



(21) 



G',G" 



3. Numerical results and discussion 



The material parameter values (p, ci, Ct) in all the materials involved (C, epoxy, 
steel, Pb and duralumin), as well as the rod/matrix parameter value contrasts 
(A/? = \p A - Pb\, Aq = \cf - cf |, Ac t = \cf - cf |), are specified in Table[T] 

The material parameter value contrasts vary substantially with structure 
composition, Ap ranging from 550 to 10157, Act from 300 to 5940, and Aq from 
672 to 10480. This variety increases the generality of our conclusions, drawn from the 
results obtained for different phononic crystals, and regarding the effect of the rod 
and lattice symmetry on the energy gap width in their elastic wave spectrum. 

In the numerical solution of lt2*U|) and lf2*T|) . 1(G) is substituted by Ic{G) f° r 
cylindrical rods (|TT))) . Irx{G) for hexagonal rods l(T7l) . Irc(G) for rectangular rods 
<|18fl and Isq(G) for square ones l|19fl . The number of reciprocal lattice vectors involved 
in Fourier expansions l|9"|l- i|ll|) is limited by condition — N < n±, n 2 < +N, confining 
the reciprocal lattice vector component values n\ and n 2 to interval [-N, +N}. From 
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Table 1. Mass density, transversal velocity and longitudinal velocity values in the 
component materials of the studied phononic crystals. Parameter value contrasts 
between the rods and the matrix (Ap, Ac; and Act) are specified as well. 





Material 


P 


Ap 


Ct 


Act 




Ac, 




[kg/m s ] 


[kg/m 3 ] 


[m/s] 


[m/s] 


[m/s] 


[m/s] 


C 


1750 




7100 




13310 




epoxy 


1200 


550 


1160 


5940 


2830 


10480 


steel 


7800 




3220 




5940 




epoxy 


1200 


6600 


1160 


2060 


2830 


3110 


Pb 


11357 




860 




2158 




epoxy 


1200 


10157 


1160 


300 


2830 


672 


duralumin 


2799 




3095 




6342 




epoxy 


1200 


1599 


1160 


1935 


2830 


3512 





Table 2. The width values of the broadest energy gap in steel/epoxy crystals 
consisting of rectangular rods of different h/h values disposed in nodes of a 
rectangular lattice (L1/L2 = 2). 



h/h 


Z mode 


XY mode 


1/2 








1 








2 


0.66 


0.50 


3 


0.47 


0.22 



the resulting finite system of algebraic equations for 2N + 1 coefficients Ui(G), 2N + 1 
allowed frequencies uj(k) are deduced for each wave vector from the first Brillouin 
zone. All the results presented in this paper were obtained for N = 7. Some dispersion 
relations were calculated for N — 10 as well, and found to be in very good conformity 
with those obtained for N — 7, regardless of the rod shape. 

Among the rod shapes considered, only the rectangular cross-section leaves some 
arbitrariness, which lies in the rectangle side ratio, Z1/Z2- In order to find the optimum 
side ratio, i.e. that corresponding to the broadest energy gap, in a rectangular 
lattice-based crystal with rectangular rods, the dependence of the gap width on the 
filling fraction was investigated for the following values of h/h- 1/2, 1, 2 and 3, 
the lattice constant ratio being fixed at Li/L 2 — 2. Some of the results obtained 
for the steel/epoxy crystal are shown in Table [3 Energy gap width values, as well 
as frequency values, are specified, throughout this study, in the units of reduced 
frequency, f2, defined as follows: 

^ = T^U (22) 

an being the distance between two neighboring lattice nodes, and c t denoting the mean 
between the transversal velocities in materials A and B. The energy gap is found to 
be the broadest for the side ratio equal to the constant lattice ratio, I1/I2 = 2 in the 
case considered. Therefore, we shall henceforth confine our analysis to rectangular 
rods with side ratio h/h — 2 and a rectangular lattice with constant lattice ratio 
L1/L2 — 2. The specific value of the constant lattice ratio, as well as that of the 
side length ratio, is of no importance for this study as long as both ratios are equal, 
which means that the symmetry of the rods corresponds to that of the lattice. Once 
the component materials, the rod shape and the crystallographic lattice type are 
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fixed, in the next step of our investigation of the effect of the symmetry relation 
between the rod and the lattice on the width of a phononic energy gap, we shall 
specify the filling fraction value corresponding to the maximum gap width. For each 
material composition, the reduced frequency is plotted against the filling fraction, /, 
on the basis of lf2T)jl and (|5TJ, with three possible lattice types (square, triangular 
and rectangular), four rod cross-section shapes (square, circle, regular hexagon and 
rectangle) and two polarizations considered; the aim is to find, in each case, the 
optimum filling fraction value, corresponding to the maximum width of the broadest 
gap. 

Figure [3 shows some examples of phononic crystal energy bands plotted against 
the filling fraction in two cases: (a) in a steel/epoxy crystal, representing the case 
of high-density rods in a low-density matrix (pg < pa), and (b) in an epoxy/steel 
structure, with low-density rods in a high-density matrix (ps > Pa)- These two 
qualitatively different cases bring out the significant effect of the material parameter 
values on the relations obtained. In the case of high-density rods in a low-density 
matrix, the optimum filling fraction value for the broadest gap (found between the 
first and the second band) is f = 0.48, the gap center shifting up the frequency scale 
as / increases (Figure EJa)). In the other case, the gaps are narrower, their centers 
shift down as / increases, and the optimum filling fraction value is much higher that 
in the case considered previously (Figure EJb)). Therefore, in the following sections 
of this paper, the cases of high-density rods in a low-density matrix and low-density 
rods in a high-density matrix will be considered separately. 

3.1. The case of high- density rods in a low-density matrix 

In this paragraph we shall focus on structures in which the matrix material is assumed 
to be epoxy; due to the lowest mass density of this material, inequality pa > Pb, 
conducive to the appearance of wider energy gaps ^SJj i s fulfilled in each case. The 
width values of the broadest energy gaps obtained in each of the examined structures, 
and corresponding to the optimum filling fraction values, are specified in Table 

The results shown in Table indicate that, among the lattice types considered 
here, it is the triangular lattice (i.e. the one whose first Brillouin zone has a shape 
closest to a circle) that allows to obtain energy gaps which are by far the broadest. 
Besides, the gaps obtained in each of the phononic crystals covered by Table E3 are 
found to be the broadest in structures in which the symmetry of the rods corresponds 
exactly to that of the lattice, which confirms the results reported in ^2]- The results 
obtained in the case of triangular lattice can be expressed by the following inequality, 
for both Z and XY modes: 



where AQ$Q is the gap width in the square-rod composite, Aflc is the corresponding 
value in the cylinder-rod structure, AQhx denotes the gap width in the composite 
with hexagonal rods and AFIrc refers to the case of rectangular rods. 

The results obtained in the case of square lattice, for both Z and XY modes, can 
be expressed as follows: 



In the case of rectangular lattice, the results obtained cannot be expressed by just one 
inequality; however, the broadest gaps are found to appear for rectangular rods. 



An HX > An c > An SQ > An RC , 



(23) 



An SQ > An c > An HX > An RC . 



(24) 
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Table 3. The width values of the broadest energy gap in the examined phononic 



crystals, in 


reduced frequency units, sq 


denotes 


square 


lattice, 


tr triangL 


rectangular 


lattice, C cylindrical 


rods, 


HX hexa 


gonal rods, SQ square i 


rectangular 


rods. 












Material 


Polarization 


Lattice 


c 


HX 


SQ 


RC 






sq 


0.37 


0.36 


0.64 


0.20 




Z 


tr 


0.77 


0.80 


0.75 


0.31 


steel 




rc 


0.04 


0.05 





0.47 


/epoxy 




sq 


0.57 


0.51 


0.74 


0.06 




XY 


tr 


0.86 


0.86 


0.75 


0.11 






rc 


0.03 








0.35 






sq 


0.21 


0.16 


0.30 


0.02 




Z 


tr 


0.49 


0.55 


0.48 


0.11 


C 




rc 





0.16 





0.27 


/epoxy 




sq 


0.13 


0.11 


0.37 







XY 


tr 


0.47 


0.49 


0.40 


o 

















0.09 






sq 


u.ou 


0.35 


0.35 


0.24 




Z 


tr 


0.47 


0.47 


0.46 


0.18 


Pb 




rc 


0.08 





0.09 


0.18 


/epoxy 




sq 


0.32 


0.32 


0.32 


0.16 




XY 


tr 


0.37 


0.37 


0.36 


0.27 






rc 











0.07 






sq 


0.21 


0.20 


0.35 


0.07 




Z 


tr 


0.48 


0.51 


0.44 


0.21 


duralumin 


rc 











0.18 


/epoxy 




sq 


0.29 


0.25 


0.37 







XY 


tr 


0.47 


0.47 


0.40 


0.06 






rc 











0.08 



The results obtained for steel/epoxy, C/epoxy and duralumin/epoxy structures 
are qualitatively the same; in the case of Z modes, all the weak inequalities in (|23|) 
and (|24(l can be replaced with strict ones. However, the results obtained for Pb/epoxy 
are different: as long as the lattice type is fixed, gaps appearing at various rod shapes 
(cylindrical, hexagonal and square) are found to be equal. The optimum filling fraction 
value lies within the range from 0.460 to 0.675 in each of the first three cases, but is 
found to be below / = 0.375, and thus relatively low, in the case of the Pb/epoxy 
composition; this means that the matrix material is then the dominant one in the 
Pb/epoxy crystal. The described properties of the energy gap width in the Pb/epoxy 
crystal are due to different contrasts of longitudinal and transversal velocity values 
in the component materials: the relations between the q and Ct values in the rods 
and in the matrix are reversed with respect to the other compositions (ctA < CtB and 
cia < cib), while the density relation is the same (pa > Pb)- The gaps are found to 
be narrowest in the C/epoxy composite, due to the lowest density contrast, Ap — 550, 
between the component materials. 

3.2. The case of low- density rods in a high-density matrix 

In order to check whether the results discussed above depend on the rod and matrix 
density values, similar computations of the phononic gap width values for different rod 
symmetries were performed for epoxy/steel, epoxy/C, epoxy/Pb and epoxy/duralumin 
crystals, representing the case of low-density rods in a high-density matrix. Providing 
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an example, Figure Hfb) shows the reduced frequency plotted against the filling 
fraction on the basis of the results obtained in the epoxy/steel crystal with cylindrical 
rods embedded in a square lattice. 

In Figure the gap width is plotted against the rod shape in steel/epoxy and 
epoxy/steel crystals (Figure (a) and (b)) as well as in Pb/epoxy and epoxy/Pb 
crystals (Figure EH (c) and (d)); Z and XY modes are considered separately (Figure 
01 (a),(c) and 01 (b),(d), respectively). The plots obtained for the remaining two 
compositions, i.e. for C/epoxy and duralumin/epoxy, are qualitatively the same as 
those obtained for the steel/epoxy crystal, and thus are not depicted in the figures 
presented in this paper. The plots shown in FigureOJindicate clearly that, regardless of 
polarization, relations i|23|) - (|24[l . obtained for crystals with high-density rods in a low- 
density matrix (represented by the solid lines), do not apply to the inverse structures 
(dashed lines), in which the density of the rod material is lower that that of the matrix. 
In triangular lattice-based crystals, hexagonal rods are found to generate narrower 
gaps than cylindrical or square rods. In square lattice-based structures, square rods 
are found to generate no gaps at all (though other rod shapes do generate them); 
also in rectangular lattice-based structures no gaps are generated by rectangular rods. 
Therefore, the inversion of the material composition clearly implies a radical change in 
the relations between the gap width and the scatterer symmetry: in none of the cases 
considered, for a given lattice type, the broadest gaps are found to appear when the 
symmetry of the rods corresponds to that of the lattice. Besides, the gaps obtained 
in such inverse structures are much narrower (four times in the epoxy/steel crystal, 
and twice in the epoxy/Pb structure) than those observed in the corresponding crystal 
with high-density rods in a low-density matrix. As similar results were obtained for all 
the compositions considered (also in the epoxy/Pb crystal, with different longitudinal 
and transversal velocity contrasts), the strongest impact on the considered relations 
can be expected to come from the mass density. 

4. Conclusion 

As indicated by the results presented in the preceding section, covering two classes 
of phononic crystals: those with high-density rods embedded in a low-density matrix, 
and those with low-density rods in a high-density matrix, material parameter values 
have much of an impact on the relation between the energy gap width and the lattice 
and rod symmetry; the effect of the mass density was found to be especially significant. 
In structures with high-density rods in a low-density matrix, the gaps are found to be 
the broadest when the symmetry of the rods corresponds to that of the lattice, i.e. 
for hexagonal rods in triangular lattice-based crystals, square rods in square lattice- 
based crystals, and rectangular rods in rectangular lattice-based crystals, as already 
demonstrated in |17] . However, this rule will cease to apply when the density of the 
rod material becomes lower than that of the matrix. In this case, when the symmetry 
of the rods corresponds to that of the lattice, gaps either fail to appear at all, or are 
much narrower than in other configurations. Gaps obtained in such inverse structures 
are always much narrower than those appearing in crystals with high-density rods in 
a low-density matrix. This is due to the fact that low-density rods embedded in a 
high-density matrix are weak scattering centers: their wave scattering effect is hardly 
noticeable and much lesser than that produced by high-density rods in a low-density 
matrix. As regards the longitudinal and transversal velocity values (q and Ct) in 
the component materials, the effect of these parameters on the investigated relations 
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proves less significant. 
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Figure captions 
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Figure 1. . Models of two-dimensional phononic crystals based on different types of 
lattice and with different rod cross-section shapes: square (sq) lattice-based crystals 
(the lattice constant being a) with rod cross section in the shape of: circle C with 
radius rrj (a), regular hexagon HX with side length b (b), square SQ with side length 
I (c), and rectangle RC with side lengths ii and I2 (d); rectangular (rc) lattice-based 
crystals (the lattice constants being L\ and L2 with rod cross section in the shape 
of: circle (e), hexagon (f), square (g) and rectangle (h); triangular (tr) lattice-based 
crystals (the lattice constant being a) with rod cross section in the shape of: circle 
(i), hexagon (j), square (k) and rectangle (I). 
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Figure 2. Z-mode reduced frequency plotted against the filling fraction, /, for a 
steel/epoxy crystal with cylindrical rods disposed in square lattice nodes (a), and its 
epoxy/steel counterpart (b). 
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Figure 3. The energy gap width (in the units of reduced frequency f2) plotted 
against the rod shape (C, HX, SQ, RC) in two-dimensional phononic crystals; (a) and 
(b) show results obtained in steel/epoxy and epoxy/steel crystals; (c) and (d) present 
those obtained for Pb/epoxy and epoxy/Pb compositions. Z and XY modes are 
considered separately, the results being depicted in (a), (c) and (b), (d), respectively. 
C denotes cylindrical rods, HX hexagonal rods, SQ square rods, and RC rectangular 
rods; the rods are disposed in nodes of a square, a triangular or a rectangular lattice, 
as indicated by symbols explained in the legend. The asterisk in the legend refers to 
the epoxy/steel structures in (a) and (b), and to epoxy/Pb structures in (c) and (d). 



